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We present a simple device based on the controlled-SWAP gate that performs quantum state 
tomography. It can also be used to determine maximum and minimum eigenvalues, expectation 
vàlues of arbitrary observables, purity estimation as well as characterizing quantum channels. The 
advantage of this scheme is that the architecture is fixed and the task performed is determined by 
the input data. 
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One of the the key issues in quantum information is, 
given an unknown quantum system, what can we learn 
about it. In particular, we are concerned not only with 
the resources needed (number of identical unknown phys- 
ical systems), but also with the complexity of quantum 
operations required (number of different devices, net- 
works, etc), in order to obtain certain information about 
a quantum state, characterized by its density matrix g. 
There are many interesting parameters of g we can deter- 
mine, such as its maximum and minimum eigenvalues, its 
purity or even g itself (state tomography [|] ) , but we also 
can use g to determine expectation vàlues of arbitrary ob- 
servables or to characterize unknown quantum channels. 
However, this usually involves building separate devices 
for each task, or even building different devices for dif- 
ferent measurements within the same task. 

In this paper we present a simple, universal device, 
whose architecture is fixed but whose behaviour is de- 
termined by the choice of input data [|| (see also || for 
a quantum optical realization of a similar idea). In fact, 
with suitable input, we can directly measure all the prop- 
erties mcntioned before. 

Consider a typical interferometric set-up for a single 
qubit: Hadamard gate, phase shift ip, Hadamard gate, 
followed by a measurement in the computational basis. 
Here and in the following, we borrow terminology from 
quantum information science and describe quantum in- 
terferometry in terms of quantum lògic gates Q. We 
modify the interferometer by inserting a controlled-Z7 op- 
eration between the Hadamard gates, with its control on 
the qubit and with U acting on a quantum system de- 
scribed by some unknown density operator p. We do 
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FIG. 1: Both the visibility and the shift of the interfer- 
ence patterns of a single qubit (top line) are affected by the 
controlled-í7 operation. 



not assume anything about the form of p, it can, for 
example, describe several entangled or separable sub- 
systems. This set-up is shown in Fig. [I]. The action of 
the controlled-í7 on p modifics the intcrference pattern 
by the factor, 



Tr pU = ve 1 ' 



(1) 



where v is the new visibility and a is the shift of the 
interference fringes, also known as the Pancharatnam 
phase [||. 

Thus, the observed visibility gives a straightforward 
way of estimating the average value of unitary operators 
U in state p and has a variety of interesting applications. 
For example, it can be used to measure some entangle- 
ment witnesses W, as long as they are unitary operators 
and the corresponding controlled-VF operations are easy 
to implement ||. Here, we focus on the applications 
related to quantum state state tomography. Clearly the 
interferometer in Fig. Q can be used to estimate any d x d 
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density matrix p by a judicious choice of d 2 — 1 unitary 
operators Í7 (the basis). However, this requires build- 
ing d 2 — 1 networks, each with a different U. We will 
show that one can build a universal state estimator i.e. a 
simple quantum device with a fixed architecture, where 
the state estimation is performed by modifying only the 
input data. 

In order to do this, let p to be composed of two sub- 
systems of dimension d, p = g a <E> Qb- We will fix our 
controllcd-ÍJ to be the controlled-y (SWAP) operator, 
such that V a) | 0) = | 0) | ot) for any pure states of the 
subsystems (Fig. g). Let us also introduce the maxi- 
mally entangled state \<fi+) = an d ^ 



|o>- H 



V 



H 



♦Pr(|o))= 



FIG. 2: Universal Quantum Estimator. By judicious choice 
of input g a , the state of Qb can be determined by examining 
the change in the visibility observed. 



.) (0+ |. In this case, the modification of the can CO nstruct the state g a 



interference pattern given by (|1|) can be written as, 
v = TV V (g a <8> Qb) = Tr g a g b . 



(2) 



Since Tr g a g b is real, we can 8x^ = and the probability 
of measuring the qubit to be in the state | 0) at the output 
is related to the visibility by, 



2Pr(| 0)) - 1. 



(3) 



If p is not separable then, writing V — P + b , we obtain 
Tr Vp = Tr p Th P + (the average value of partially trans- 
posed p in the maximally entangled state P+). 

Now, let g b be an unknown dxd density operator. Such 
an operator is determined by d 2 — 1 real parameters. In 
order to estimate matrix elements of g b in a prescribed 
orthonormal basis, { | n) } , we proceed as follows: We run 
the interferometer as many times as possible (limited by 
the number of copies of g b at our disposal) on the input 
\ip) {ip | ® Qb, where U/>) is a pure state of our choice. For 
a fixed | tp), after several runs we obtain an estimation 
of, 



V= Q b \ tp). 



(4) 



The diagonal elements (n | g b | n) can be determined us- 
ing the input states | ra) {n | (8 Qb- The real part of the 
off-diagonal element (n \ g b \ k) can be estimated by choos- 
ing | tfj) = (I n) + I fe))/\/2, and the imaginary part by 
choosing \ ip) = (| n) + i\ fc))/v / 2· In particular, if we 
want to estimate the density operator of a qubit, we can 
choose the pure states, | 0) (spin +z), (| 0) + | l)) /\/2 
(spin +x) and (| 0) + i\ l)) /y/2 (spin +y). 

We can extend the procedure above to estimate expec- 
tation vàlues of arbitrary observables. It can be shown 
that the mean value of an arbitrary observable can be 
reduced to the estimation of a binary two-output POVM 
0. Similarly, the mean value (A) 8b of an arbitrary ob- 
servable A in state g b can be measured using the setup 
in Fig. ^| with a suitable input g a ■ We shall apply the 
technique utilized in Refs. [|[ B. As A' — 7I + A is pos- 
itive if —7 is the minimum negative eigenvalue of A, we 



QA ' = Tr (A>) and apply our 
interference scheme to the pair g^> ® g b . The visibility 
gives us the mean value of V (SWAP), 

v = ( v )g A ^et = Tr(gA>gb), (5) 

which lcads us to the desired value, 

(A) el· =Tr(Q b A)=vTrA + 1 (vd-l), (6) 

where Tr I = d. 

Other quantities related to g b may be determined by 
simple modification of the interferometry scheme. If we 
have at our disposal two copies of gb per run, by running 
the interferometer on the input, p = gb®g bl the resulting, 



(7) 



where {A^} are the eigenvalues of g b , gives us an estimate 
of the purity of g b . In the single qubit case, this mea- 
surement allows us to estimate various functionals of g b , 
such as the length of the Bloch vector, 
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(8) 



Note that the direction, however, is left completely 
undetermined. The procedure of estimating eigenval- 
ues and non-linear functionals of g b can be generalized 
for larger dimensional systems, but requires controlled- 
SHIFT gates g. 

By adapting the input to the interferometer, we can 
also estimate the extremal eigenvalues and eigenvectors 
of g b . In this case, the input states are also of the form 
\tl>) (ip \ ®Q b but we vary | ip) and search for the minimum 
and the maximum of v = {tp | Qb \ ip)- This is usually a 
complicated task as it involves scanning 2(d~ 1) param- 
eters of tp. The visibility is related to the overlap of the 
reference state, | tp) and gb by, 

V4, = Tr ^| Tp) (ip | Ç Ai | r/,) (77, |J 

= E A *K^>i 2 = E A ^' ( 9 ) 
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where YliPi — !• This is a convex sum of the eigenval- 
ues of Qb and is minimized (maximized) when | -0) = 
| r)min) (\vmax))- For any \ ip) ^ \i] mul ) (\r) max )), 
there exists a state, | -0'), in the neighbourhood of | ip) 
such that v^j/ < (v^/ > v^,) thus this global opti- 
mization problem can be solved using Standard iterative 
methods, such as steepest decent JÏ0[ . 

Estimation of extremal eigenvalues plays a significant 
role in the direct detection of quantum entanglement ||] 
and distillation |TÏ| . In some special cases of two qubits 
described by the density operator gi, such that at least 
one of the qubits is in the maximally mixed state, we 
can test for the separability of gi, by checking whether 
the maximal eigenvalue of gb does not exceed \ (we shall 
return to this special case in our subsequent discussion 
of the quantum channel tomography). 

Let us now turn our attention to characterizing quan- 
tum channels. Recali that a quantum channel is a trace 
preserving linear map, g — > A(g), which takes quantum 
states to quantum states, and whose trivial extensions, 
Ifc ® A do the same, i.e. A is a complctcly positive map. 
Using the well known Jamiolkowski isomorphism fÏ2f be- 
tween quantum channels and bipartite states, quantum 
channels can be characterized in a simple way The single 
qubit channel capacity and the distillability of entangled 
states can be determined by extremal eigenvalue estima- 
tion. 

Suppose we have an unknown quantum channel, A, 
which we would like to characterize. The Jamiolkowski 
isomorphism identifies a quantum channel with its action 
on half of a maximally entangled state. The procedure is 
thus (Fig. ||): We prepare maximally entangled states of 
two partides P+ = \ J2ij | *) I ® | *) C? |j We send one 



Particle 1 



particle through the channel, 

P + -»• [I ® A] P + 



pa; 



(10) 



We then estimate 



which now characterizes A. We intèrpret this as A map- 
ping the | i) (j | í,l -element of an input density matrix to 
the output matrix, A (I i) (j |). Thus, knowledge of g\ 
allows us to determine the action of A on an arbitrary 
state, g — > A(g). 



Now consider the case where A is a single qubit chan- 
nel. The channel capacity, Q(A) is an important param- 
eter of A [|Ï3| This is the optimal rate of reliably 

sending qubits per use of A. Depending on whether Al- 
ice and Bob can send classical information to each other 
as an additional resource, one can consider several capac- 
ities, Qc where C = 0, <—,—>, <->, corresponding to zero 
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FIG. 3: A quantum channel is isomorphic with its action on 
half of a maximally entangled bipartite state. 



way, one way and two way classical communication. In 
general, it is very difficult to calculate the capacity of a 
given channel, A. Here we shall provide a simple neces- 
sary and sufficient condition for a one qubit channel to 
have non-zero two-way capacity, > (obviously a 
necessary condition for the other three capacities to be 
non-zero), which can be determined using our extremal 
eigenvalue estimation scheme. 

In order to do this, we shall use the state, qa, defined 
111 ©. Clearly, for any channel A, gA is maximally 
mixed when reduced to the subsystem A. It is also true 
that for any state gA with a maximally mixed subsystem 
A there exists a channel A which generates gA via the 
formula (|ï^). If gA is maximally mixed when reduced to 
both subsystems A and B, then the channel A is called 
bistochastic - it maps maximally mixed states into max- 
imally mixed states. (Two-qubit states, corresponding 
to bistochastic channels, have been completely charac- 
terized geometrically as tetrahedrons in R 3 |Ï5|, |Ï6|]). 

It is known (sec |Ï7| ) that a two qubit state is two- 
way distillable iff the operator gA® I — gAB has negative 
eigenvalue. Now for states of the type gA (those are all 
states with gA = 5), this reduces to the requirement 
that gA has maximal eigenvalue greater than i. This is 
also equivalent to Q W (A) > 0, since two-way distillable 
entanglement (which is non-zero iff given state is two way 
distillable) is simply the lower bound for Q^(A) 14 . 



Now we can apply our estimator of the maximal eigen- 
value to check whether a given one-qubit channel has a 
non-zero two-way capacity. Instead of fully determining 
the channel as in the previous section, we simply estimate 
the maximum eigenvalue of gA- This scheme can be also 
be used to find whether a given two qubit state, with one 
sub-system maximally entangled, is two-way distillable. 

We have presented a universal quantum estimator, 
whose action is determined by the input. The is based 
on interferometry and the controlled-SWAP operation. 
By suitable choice of input states, we are able to per- 
forin quantum tomography, extremal eigenvalue estima- 
tion, purity tests and quantum channel characterization. 
Finally let us mention that the controlled-SWAP opera- 
tion is a direct generalization of a Fredkin gate pll and 
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can be constructed out of simple quantum lògic gates [|Ï9| . 
This means that experimental realizations of a universal 
quantum estimator are within thc rcach of quantum tech- 
nology that is currently being developed. 
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